Análisis de la variabilidad genética de una muestra de la población de Bogotá: hacia la constitución de un mapa de haplotipos / Genetic variability analysis in a population from Bogotá: Towards a haplotype map

Resumen Introducción. Los proyectos del mapa de haplotipos (HapMap) y de los 1.000 genomas han sido fundamentales para la compresión del componente genético de las enfermedades comunes y los fenotipos normales. Sin embargo, la variabilidad genética colombiana incluida en estos proyectos no es representativa del país. Objetivo. Contribuir al conocimiento de la variabilidad genética de la población colombiana a partir del estudio genómico de una muestra de individuos de Bogotá. Materiales y métodos. Se genotipificaron 2'372.784 marcadores genéticos de 32 individuos nacidos en Bogotá y de padres originarios de la misma ciudad utilizando la plataforma Illumina™. Los niveles de variabilidad genética se determinaron y se compararon con los datos disponibles de otras poblaciones del proyecto de los 1.000 genomas. Resultados. Los individuos analizados presentaron una variabilidad genética semejante a la de poblaciones con las que comparten ancestros. No obstante, a pesar de la poca diferenciación genética detectada en la población de Bogotá y en la de Medellín, el análisis de los componentes principales sugiere una composición genética diferente en las dos poblaciones. Conclusiones. El análisis genómico de la muestra de Bogotá permitió detectar similitudes y diferencias con otras poblaciones americanas. El aumento de tamaño de la muestra bogotana y la inclusión de muestras de otras regiones del país permitirán una mejor compresión de la variabilidad genética en Colombia, lo cual es fundamental para los estudios de salud humana, y la prevención y el tratamiento de enfermedades comunes en el país.

Abstract Introduction: The HapMap and the 1000 Genomes projects have been important for understanding the genetic component of common diseases and normal phenotypes. However, the Colombian genetic variability included in these projects is not fully representative of our country. Objective: To contribute to the knowledge of the Colombian genetic variability through the genomic study of a sample of individuals from Bogotá. Materials and methods: A total of 2,372,784 genetic markers were genotyped in 32 individuals born in Bogotá whose parents are from the same region, using the Illumina™ platform. The genetic variability levels were determined and compared with the data available from other populations of the 1000 Genomes Project. Results: The genetic variability detected in the individuals from Bogotá was similar to those with shared ancestry. However, despite the low levels of genetic differentiation between Bogotá and Medellín, populations the principal component analysis suggested a different genetic composition in them. Conclusions: Our genomic analysis of a Bogotá sample allowed us to detect similarities and differences with other American populations. The increase of the Bogotá sample and the inclusion of samples from other regions of the country will improve our understanding of the genetic variability in Colombia, essential for studies of human health and the prevention and treatment of common diseases in our country.

Semantic Modeling for SNPs Associated with Ethnic Disparities in HapMap Samples

Single-nucleotide polymorphisms (SNPs) have been emerging out of the efforts to research human diseases and ethnic disparities. A semantic network is needed for in-depth understanding of the impacts of SNPs, because phenotypes are modulated by complex networks, including biochemical and physiological pathways. We identified ethnicity-specific SNPs by eliminating overlapped SNPs from HapMap samples, and the ethnicity-specific SNPs were mapped to the UCSC RefGene lists. Ethnicity-specific genes were identified as follows: 22 genes in the USA (CEU) individuals, 25 genes in the Japanese (JPT) individuals, and 332 genes in the African (YRI) individuals. To analyze the biologically functional implications for ethnicity-specific SNPs, we focused on constructing a semantic network model. Entities for the network represented by "Gene," "Pathway," "Disease," "Chemical," "Drug," "ClinicalTrials," "SNP," and relationships between entity-entity were obtained through curation. Our semantic modeling for ethnicity-specific SNPs showed interesting results in the three categories, including three diseases ("AIDS-associated nephropathy," "Hypertension," and "Pelvic infection"), one drug ("Methylphenidate"), and five pathways ("Hemostasis," "Systemic lupus erythematosus," "Prostate cancer," "Hepatitis C virus," and "Rheumatoid arthritis"). We found ethnicity-specific genes using the semantic modeling, and the majority of our findings was consistent with the previous studies - that an understanding of genetic variability explained ethnicity-specific disparities.

A Short History of the Genome-Wide Association Study: Where We Were and Where We Are Going

Recent rapid advances in genetic research are ushering us into the genome sequence era, where an individual's genome information is utilized for clinical practice. The most spectacular results of the human genome study have been provided by genome-wide association studies (GWASs). This is a review of the history of GWASs as related to my work. Further efforts are necessary to make full use of its potential power to medicine.

Population analysis of vitamin D receptor polymorphisms and the role of genetic ancestry in an admixed population

The vitamin D receptor (VDR) is an essential protein related to bone metabolism. Some VDR alleles are differentially distributed among ethnic populations and display variable patterns of linkage disequilibrium (LD). In this study, 200 unrelated Brazilians were genotyped using 21 VDR single nucleotide polymorphisms (SNPs) and 28 ancestry informative markers. The patterns of LD and haplotype distribution were compared among Brazilian and the HapMap populations of African (YRI), European (CEU) and Asian (JPT+CHB) origins. Conditional regression and haplotype-specific analysis were performed using estimates of individual genetic ancestry in Brazilians as a quantitative trait. Similar patterns of LD were observed in the 5' and 3' gene regions. However, the frequency distribution of haplotype blocks varied among populations. Conditional regression analysis identified haplotypes associated with European and Amerindian ancestry, but not with the proportion of African ancestry. Individual ancestry estimates were associated with VDR haplotypes. These findings reinforce the need to correct for population stratification when performing genetic association studies in admixed populations.

A Scheme for Filtering SNPs Imputed in 8,842 Korean Individuals Based on the International HapMap Project Data

Genome-wide association (GWA) studies may benefit from the inclusion of imputed SNPs into their dataset. Due to its predictive nature, the imputation process is typically not perfect. Thus, it would be desirable to develop a scheme for filtering out the imputed SNPs by maximizing the concordance with the observed genotypes. We report such a scheme, which is based on the combination of several parameters that are calculated by PLINK, a popular GWA analysis software program. We imputed the genotypes of 8,842 Korean individuals, based on approximately 2 million SNP genotypes of the CHB+JPT panel in the International HapMap Project Phase II data, complementing the 352k SNPs in the original Affymetrix 5.0 dataset. A total of 333,418 SNPs were found in both datasets, with a median concordance rate of 98.7%. The concordance rates were calculated at different ranges of parameters, such as the number of proxy SNPs (NPRX), the fraction of successfully imputed individuals (IMPUTED), and the information content (INFO). The poor concordance that was observed at the lower values of the parameters allowed us to develop an optimal combination of the cutoffs (IMPUTED> or =0.9 and INFO> or =0.9). A total of 1,026,596 SNPs passed the cutoff, of which 94,364 were found in both datasets and had 99.4% median concordance. This study illustrates a conservative scheme for filtering imputed SNPs that would be useful in GWA studies

Chromosome 22 LD Map Comparison between Korean and Other Populations

Single nucleotide polymorphisms (SNPs) are the most abundant forms of human genetic variations and resources for mapping complex genetic traits and disease association studies. We have constructed a linkage disequilibrium(LD) map of chromosome 22 in Korean samples and compared it with those of other populations, including Yorubans in Ibadan, Nigeria (YRI), Centred'Etude du Polymorphisme Humain (CEPH) reference families (CEU), Japanese in Tokyo (JPT) and Han Chinese in Beijing (CHB) in the HapMap database. We genotyped 4681 of 111,448 publicly available SNPs in 90 unrelated Koreans. Among genotyped SNPs, 4167 were polymorphic. Three hundred and five LD blocks were constructed to make up 18.6% (6.4 of 34.5 Mb) of chromosome 22 with 757 tagSNPs and 815 haplotypes(frequency > or = 5.0%). Of 3430 common SNPs genotyped in all five populations, 514 were monomorphic in Koreans. The CHB + JPT samples have more than a 72% overlap with the monomorphic SNPs in Koreans, while the CEU + YRI samples have less than a 38% overlap. The patterns of hot spots and LD blocks were dispersed throughout chromosome 22, with some common blocks among populations, highly concordant between the three Asian samples. Analysis of the distribution of chimpanzee-derived allele frequency (DAF), a measure of genetic differentiation, Fst levels, and allele frequency difference (AFD) among Koreans and the HapMap samples showed a strong correlation between the Asians, while the CEU and YRI samples showed a very weak correlation with Korean samples. Relative distance as a quantitative measurement based upon DAF, Fst, and AFD indicated that all three Asian samples are very proximate, while CEU and YRI are significantly remote from the Asian samples. Comparative genome-wide LD studies provide useful information on the association studies of complex diseases.

Chromosome 22 LD Map Comparison between Korean and Other Populations

Single nucleotide polymorphisms (SNPs) are the most abundant forms of human genetic variations and resources for mapping complex genetic traits and disease association studies. We have constructed a linkage disequilibrium(LD) map of chromosome 22 in Korean samples and compared it with those of other populations, including Yorubans in Ibadan, Nigeria (YRI), Centred'Etude du Polymorphisme Humain (CEPH) reference families (CEU), Japanese in Tokyo (JPT) and Han Chinese in Beijing (CHB) in the HapMap database. We genotyped 4681 of 111,448 publicly available SNPs in 90 unrelated Koreans. Among genotyped SNPs, 4167 were polymorphic. Three hundred and five LD blocks were constructed to make up 18.6% (6.4 of 34.5 Mb) of chromosome 22 with 757 tagSNPs and 815 haplotypes(frequency > or = 5.0%). Of 3430 common SNPs genotyped in all five populations, 514 were monomorphic in Koreans. The CHB + JPT samples have more than a 72% overlap with the monomorphic SNPs in Koreans, while the CEU + YRI samples have less than a 38% overlap. The patterns of hot spots and LD blocks were dispersed throughout chromosome 22, with some common blocks among populations, highly concordant between the three Asian samples. Analysis of the distribution of chimpanzee-derived allele frequency (DAF), a measure of genetic differentiation, Fst levels, and allele frequency difference (AFD) among Koreans and the HapMap samples showed a strong correlation between the Asians, while the CEU and YRI samples showed a very weak correlation with Korean samples. Relative distance as a quantitative measurement based upon DAF, Fst, and AFD indicated that all three Asian samples are very proximate, while CEU and YRI are significantly remote from the Asian samples. Comparative genome-wide LD studies provide useful information on the association studies of complex diseases.

Development of KHapmap Browser using DAS for Korean HapMap Research

The Korean HapMap Project has been carried out for the last 5 years since it started in June, 2003. The project generated data for a sum of 1,764,000 Korean SNPs and formally registered the data to the dbSNP of NCBI (The dbSNP website. 2008). We have developed a series of software programs for association studies as well as for the comparison and analysis of Korean HapMap data with four other populations (CEPH, Yoruba, Han Chinese, and Japanese populations). The KHapmap Browser was developed and integrated to provide haplotype retrieval and comparative study tools of human ethnicities for comprehensive disease association studies (http://www.khapmap.org). On that basis, GBrowse was adopted in the KHapmap Browser for inherent Korean genetic data, and a provision of extended services was pledged with the distributed sequence annotation system (DAS). The dynamic linking service of the KHapmap Browser to other tools in our intranetwork environment provides many enhanced functions over GBrowse without DAS. KHapmap Browser is expected to be an invaluable tool for the study of Korean and international Hapmap data.

The Korean HapMap Project Website

Single nucleotide polymorphisms (SNPs) are the most abundant form of human genetic variation and are a resource for mapping complex genetic traits. A genome is covered by millions of these markers, and researchers are able to compare which SNPs predominate in people who have a certain disease. The International HapMap Project, launched in October, 2002, motivated us to start the Korean HapMap Project in order to support Korean HapMap infrastructure development and to accelerate the finding of genes that affect health, disease, and individual responses to medications and environmental factors. A Korean SNP and haplotype database system was developed through the Korean HapMap Project to provide Korean researchers with useful data- mining information about disease-associated biomarkers for studies on complex diseases, such as diabetes, cancer, and stroke. Also, we have developed a series of software programs for association studies as well as the comparison and analysis of Korean HapMap data with other populations, such as European, Chinese, Japanese, and African populations. The developed software includes HapMapSNPAnalyzer, SNPflank, HWE Test, FESD, D2GSNP, SNP@Domain, KMSD, KFOD, KFRG, and SNP@WEB. We developed a disease-related SNP retrieval system, in which OMIM, GeneCards, and MeSH information were integrated and analyzed for medical research scientists. The kHapMap Browser system that we developed and integrated provides haplotype retrieval and comparative study tools of human ethnicities for comprehensive disease association studies (http://www.khapmap.org). It is expected that researchers may be able to retrieve useful information from the kHapMap Browser to find useful biomarkers and genes in complex disease association studies and use these biomarkers and genes to study and develop new drugs for personalized medicine.

Haplotype Phylogeny of a 200kb Region in the Human Chromosome X Terminal Band (q28)

The haplotypes of a 200 kb region in the human chromosome X terminal band (q28) were analyzed using the International HapMap Project PhaseII data, which had been collected for three analysis panels (YRI, CEU, and CHB+JPT). When multiple linkage disequilibrium blocks were encountered for a panel, the neighboring haplotypes that had crossover rate of 5% or more in the panel were combined to generate 'haploid' configurations. This resulted in 8, 7, and 5 'haploid' configurations for the panels of YRI, CEU, and CHB+JPT, respectively. The multiple sequence alignment of these 'haploids' was used for the calculation of allele-sharing distances and the subsequent principal coordinate analysis. Two 'haploids' in CEU and CHB+JPT were hypothesized as 'parental' in light of the observations that the successive recombinants of these haploids can model two other haploids in CEU and CHB+JPT, and that their configurations were consistent with those in YRI. This study demonstrates the utility of haplotype phylogeny in understanding population evolution.

FESD II: A Revised Functional Element SNP Database of Human Ethnicities

The Functional Element SNPs Database (FESD) categorizes functional elements in human genic regions and provides a set of single nucleotide polymorphisms (SNPs) located within each area. Users may select a set of SNPs in specific functional elements with haplotype information and obtain flanking sequences for genotyping. Our previous version of FESD has been improved in several ways. We regenerated all the data in FESD II from recently updated source data such as HapMap, UCSC GoldenPath, dbSNP, OMIM, and TRANSFAC(R). Users can obtain information about tagSNPs and simulate LD blocks for each gene from four ethnicities in the HapMap project on the fly. FESD II employs a Java/JSP web interface for better platform portability and higher speed than PHP in the previous version. As a result, FESD II provides its users with more powerful information about functional element SNPs of human ethnicities.